Reentrant stability of superconducting films 
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We propose a mechanism responsible for the abrupt vanishing of the dendritic flux instability found 
in many superconducting films when an increasing magnetic field is applied. The onset of flux 
avalanches and the subsequent reentrance of stability in NbN films was investigated using magneto- 
optical imaging, and the threshold fields were measured as functions of critical current density, 
j c - The results are explained with excellent quantitative agreement by a thermomagnetic model 
published recently, Phys. Rev. B 73, 0145f2 (2006), showing that the reentrant stability is a direct 
consequence of a monotonously decreasing j c versus field. 
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In most superconductors a slow increase of external 
magnetic field is accompanied by a gradual penetration 
of magnetic vortices, which make their way into the sam- 
ple through a random energy landscape created by struc- 
tural irregularities. This results in formation of a critical 
state characterized by a critical gradient of the vortex 
density corresponding to the maximum lossless current 
in the superconductor. [l[ The critical state, however, is 
metastable and can be destroyed by flux jumps where 
large-scale redistribution of the vortices suddenly takes 
place. Such dramatic events, observed experimentally as 
abrupt drops in the magnetization, are due to a thermo- 
magnetic instability where the local dissipation associ- 
ated with vortex motion reduces the pinning, which in 
turn facilitate further motion. Q With this positive feed- 
back a small perturbation can quickly evolve into a flux 
avalanche of sample spanning dimensions. 

In recent years, space-resolved measurements, in par- 
ticular using magneto-optical imaging (MOI), have re- 
vealed the detailed morphology of flux avalanches. Most 
work was carried out on superconducting thin films of 
MgB 2 J,|,[iJi0L Nb,@,i Pb,0 Nb 3 Sn,dH NbN, [3 
YBa 2 Cu 3 O x ,(ll[l3 and YNi 2 B 2 C,[l5l where the images 
show that a typical flux avalanche in superconducting 
films has a branched dendritic structure with 50-100 mi- 
cron wide branches. While each avalanche event has 
only local impact, they occur very frequently, i.e., with 
very small applied field intervals, and therefore manifest 
in magnetization versus field curves as a strong noisy 
component. [1 S, E3, El H E3] 

These magnetization curves display also two other 
generic features, namely the existence of a threshold field 
for the onset of avalanche activity, typically at a few mil- 
litesla, and an upper threshold field above which the su- 
perconductor regains full stability. In a recent Letter 
we reported a detailed study of the lower threshold, and 



explained key features of the instability onset. The im- 
portant question why superconducting films suddenly be- 
come stable again at high fields remains open, and is the 
focus of the present work. We report here experiments 
carried out on films of superconducting NbN, and discuss 
the results within the theoretical framework of a recently 
proposed thermo-magnetic model H 2^] ■ The reentrant 
stability is explained by accounting for the field depen- 
dence of the critical current density. 

Superconducting films of NbN were prepared on single 
crystal (100) MgO substrates by pulsed laser ablation of 
a high purity Nb target in a controlled ultra high purity 
(99.9996) N 2 environment. A KrF excimer laser (248 nm) 
operated at 20 Hz was used for ablation with pulse energy 
density of 5 J/cm 2 on the target surface. Detailed pro- 
cedures and structural characterization of the samples, 
which are superconducting below T c = 15 K, is found in 
Ref.~ 
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The films have a thickness of d = 280 nm, and 
were lithographically patterned into a rectangular shape 
2.4 mm wide and 4.8 mm long using reactive ion etching 
in a CHF3 + 2 plasma. 

Flux penetration into the samples was investigated us- 
ing a MOI method based on the Faraday effect in an 
in-plane magnetized bismuth-substituted ferrite garnet 
sensor film, that was placed on top of the NbN film. The 
sample was mounted on the cold finger of a continuous 
He-flow cryostat (Microstat, Oxford), which in turn was 
placed under the objective lens of a reflected light po- 
larizing microscope (Leica). A variable magnetic field 
was applied perpendicular to the sample using a pair of 
resistive coils. The microscope was equipped with a mer- 
cury lamp, and by using crossed polarizers the resulting 
image brightness represents the local magnitude of the 
perpendicular flux density. 

Shown in Fig. [T] is a set of images recorded during a 
slow ramping up of the applied magnetic field, H, after 
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FIG. 1: Magneto-optical images of a NbN film (only left half 
of the sample area is shown) at 4 K in increasing applied 
field. The left edge of the sample, seen as a vertical line 
of enhanced brightness, is 2.4 mm long. The images were 
contrast enhanced individually. 



the sample was initially zero-field cooled to T = 4 K. In 
image (a), taken at H = 8 Oe, there is only a shallow pen- 
etration, with a flux front that moved gradually inwards 
as the field increased. This smooth mode of penetration 
ended abruptly when reaching H = 12 Oe, where a flux 
avalanche suddenly occurred. The avalanche area is seen 
in Fig.Q]taken at 12 Oe as the branching structure rooted 
on the lower sample edge. As the applied field increased 
further, more and more avalanches took place, creating 
a complex pattern of flux dendrites covering most of the 
sample area, see Fig. Q] at H = 42.5 Oe. Then, as the 
field reached H = 357 Oe, the avalanche activity stopped 
entirely although the field continued to increase. The ad- 
vancing flux front now erased the previously formed den- 
dritic structures, see Fig. Q] at H = 400 Oe and 425 Oe, 
and the critical state mode of flux penetration was re- 
stored. The bright fan-like features in the lower part of 
the images are due to tiny film defects, and are not re- 
lated to the flux instability. 




Critical current density, j c 



FIG. 2: Schematic plot illustrating the existence of two 
threshold fields for the dendritic instability. The two main 
curves represent the threshold field H th versus j c according 
to Eqs. ([T]) and ([2]), and a typical monotonous field depen- 
dence of the critical current density. The intersection of the 
two curves defines the two thresholds H{ h and H2 1 for the 
onset and vanishing of the instability. 

This experiment gives direct evidence that there ex- 
ists an upper threshold field for the dendritic instabil- 
ity, in agreement with earlier measurements on various 
superconducting films. 0, [H, [H, [lj| Before presenting 
more results, we give first a qualitative argument for why 
reentrant stability can follow from the thermo-magnctic 
model proposed in Ref. [2(1 The theory predicts that a 
superconducting thin strip of half-width w becomes un- 
stable when the flux penetration depth exceeds 

„„ 7T IkT* ( hh Q T*\ 

provided that I* < w. Here, j c is the critical current 
density, T* = -(din jc/dT)" 1 , E is the electric field, 
k is the thermal conductivity, and ho is the coefficient 
of heat transfer from the superconducting film to the 
substrate. The parameter n characterizes the nonlincar- 
ity of the current-voltage curve of the superconductor, 
n = din E/ din j ^S> 1. The instability onset field can be 
obtained from £* using the Bean model relation between 
the penetration depth and applied field. [2rj| 

#th = 3cd gxfxoA (-^—) . ( 2 ) 

7T \W — I* J 

It follows from Eqs. {T]) and |J2J) that the threshold field 
depends strongly on the critical current density, more 
specifically according to the graph marked H th (j c ) in 
Fig. [2] For intermediate j c the threshold field is nearly 
constant, and increases slowly as j c becomes larger giv- 
ing eventually a sublinear asymptotic dependence, H th ~ 
jc^ 4 ■ More importantly, when j c decreases the threshold 
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FIG. 3: Field dependence of the critical current density at 
T = 4K measured using MOI. 



FIG. 4: The lower and upper threshold fields Hf 1 (disks) 
and H2 (triangles) measured at different frozen-in fields after 
field-cooling to 4 K. 



field will diverge at some finite j c , corresponding to (*{j c ) 
approaching w. Taking now into account that the critical 
current density most commonly decreases with the field, 
typically as the j c {B) curve in Fig. [2j it follows that the 
field range with unstable behavior can indeed have both 
a lower and upper limit, iJj h and H^ 1 , as indicated in 
the figure. Evidently, it is essential in this picture how 
fast j c is decreasing with the field. E.g., in the Bean 
model (constant j c ), or with a j c having only weak field 
dependence, the conditions for having an upper thresh- 
old will never be met. Note also that if j c is sufficiently 
small, e.g., because of a larger temperature, the thcrmo- 
magnetic avalanches will not occur at any magnetic field. 

To obtain quantitative support for this explanation 
we performed MOI of flux penetration after first cooling 
the sample to 4 K in various constant magnetic fields, 
Hf c . This allowed measurement of the field dependence 
of both j c and the two instability threshold fields. Cool- 
ing fields up to 300 Oe were applied, and MOI did not 
detect any contrast, which implies that the full flux den- 
sity H{ c was frozen into the sample. To measure j c a 
small additional field was subsequently applied, crea ting 
a critical state type of penetration of the new flux. [27J 
From the depth of the penetration front in the middle sec- 
tion of the rectangular sample the j c was determined, [26| 
and Fig. [3] shows its dependence on the cooling field. A 
nearly exponential decay was found, with j c decreasing 
from 7.6 • 10 10 Am -2 at zero field to almost one half at 
Hfc = 300 Oe. This agrees well with the results obtained 
earlier from ac susceptibility measurements [2ll ] on the 
same type of films. 

The threshold fields were determined by slowly ramp- 
ing up the additional field after an initial field-cooling. 
The values of the total applied field when the first and 
last flux dendrite appeared were recorded for various Hf c , 
giving the result shown in Fig. 2] The onset field in- 
creases almost linearly, and with some upward curvature, 



whereas the upper threshold field remains essentially con- 
stant until the two thresholds eventually merge into one 
when Hf c = 300 Oe. For larger cooling fields, i.e., for 
smaller j c 's, the instability was found to be fully sup- 
pressed. 

To demonstrate that these data are in excellent agree- 
ment with the scenario depicted in Fig. [2J we note 
first that the applied field which induces shielding cur- 
rents in the sample, and triggers the first avalanche, is 
AH = H[ h — Hf c . Thus, in the presence of a cool- 
ing field, the curve representing the H th (j c ) in Fig. [2] 
should be shifted upwards by the amount Hf c , while the 
jc(B) curve remains as is. Then, as Hf c becomes larger 
the distance between the two intersection points of the 
curves becomes gradually smaller, and at some cooling 
field the two thresholds become one, just like in the ex- 
periments. To show the agreement quantitatively, the 
data for the lower threshold field are replotted in Fig. 
as AH versus j c - The full curve is a theoretical curve 
obtained from Eqs. (p} and using two fitting pa- 
rameters, kT*/E = 343.75 W/V and h a T*/nE = 687.5 
kW/V (which can mean, e.g., n = 5 W/Km, ho = 110 
kW/Km, n = U,E= 0.16 V/m, T* = 11 K ). 

In contrast to H\ , the upper threshold field flj was 
found to be essentially independent of Hf c . Even this 
behavior is fully consistent with our explanation illus- 
trated by Fig. [21 Near the intersection point determining 
H2 1 , the diverging H th (j c ) curve has a very steep slope. 
Hence, a vertical shift of H th (j c ) accounting for a non- 
zero cooling field Hf c hardly changes the value of H\ . 
The robustness of the upper threshold field by being in- 
dependent of the initial conditions, is actually a remark- 
able and unexpected result. Indeed, the flux distribution 
after cooling the sample in Hf c and then ramping the field 
up to H% dramatically depends on the value of Hf c . In 
the zero- field cooled experiment shown in Fig. [T] one finds 
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FIG. 5: Instability onset field AH as a function of the critical 
current at 4 K. The plot shows experimental data (squares) 
and a theoretical prediction (full line) using fitting parameters 
given in the text. This plot made use also of the experimental 
data in Fig. [3] to convert from Hf c to j c . The error bars 
indicate the spread in the data obtained by repeating the 
experiments several times. 



dozens of overlapping dendrites, while for Hf c = 300 Oe 
the flux distribution remains of a critical-state type un- 
til a single dendritic structure is formed at flj • Note 
that the exact dendritic pattern is every time different 
even if one repeats the experiment under the same con- 
ditions. Thus, despite the wide variety of flux distribu- 
tions, the dendritic activity in all these cases vanishes at 
the same threshold field Hf w 350 Oe. Note also that 
the H th (j c ) curve in Fig. [2] was calculated [2oT] assuming 
an initial critical-state flux distribution, and therefore 
the observed invariance of H^ 1 suggests that the theoret- 
ical prediction is valid for a much broader set of initial 
magnetic states. 

The present thermomagnetic scenario, and in partic- 
ular the H th (j c ) dependence shown in Fig. has very 
recently received even further support from MOI studies 
of MgB2 films grown on vicinal substrates giving a slight 
anisotropy in the critical current. [25| The observed dra- 
matic anisotropy in the direction of dendritic avalanches 
was fully explained by the very steep slope of the H th (j c ) 
curve. 

In conclusion, the dendritic flux instability, which is a 
serious threat to superconducting devices based on thin 
films, was investigated to understand mechanisms that 
allow the films to recover stability. We explained reen- 
trant stability in the framework of a thermomagnetic 
model, [2(| and demonstrated by a quantitative agree- 
ment of the theory with experimental results, that it 
stems from the field dependence of the critical current. 
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